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■ relativistic case. Moreover, these solutions are uniquely characterised as fixed points of 

a certain operator. We also find a global energy functional for the system allowing us to 
obtain L p -nonlinear stability of some particular equilibria in the periodic setting. 
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1. Introduction 

Given a population of electrons, with mass m and charge — e, assumed to be rela- 
tivistic, we denote 

P P 



We analyse a reduced ID Vlasov— Maxwell system introduced recently in the physical 
literature for studying laser-plasma interaction. This system can be seen as a standard 
Vlasov equation in which the field is split in two terms: an electrostatic field obtained 
from Poisson's equation and a vector potential term satisfying a nonlinear wave equation. 
Both nonlinearities in the Poisson and wave equations are due to the coupling with the 
Vlasov equation through the charge density. We show global existence of weak solutions 
in the non-relativistic case, and global existence of characteristic solutions in the quasi- 
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the velocity corresponding to a given momentum p. The electrons move under the 
effect of an electric field E and a magnetic field B. Then, their distribution function 
f(t, q, p), where q denotes the position variable, is solution to the Vlasov equation: 

| + v(p).g-e(E + v(p)xB).|Uo. (1.1) 
The fields E and B are the sum of three parts: 

(1) the self-consistent fields created by the electrons; 

(2) the electromagnetic field of a laser wave which is sent into the medium (called 
the pump wave); 

(3) the electrostatic field E ex t(q) generated by a background of ions which are 
considered immobile during the time scale of the wave, and/or by an external, 
static confinement potential. 

In all cases, we denote by n cxt := £o divE 0Xt /e. Without this term, the population 
of electrons could not be dynamically stable. Then, the Maxwell system is written: 

f.^B + ij, (1.2) 

— = -curlE, (1.3) 
divE = — (riext-n), (1.4) 
divB = 0, (1.5) 

where c and So are the speed of light and the dielectric permittivity of vacuum, and 
the electron density and flux n and j are the first two moments of the distribution 
function f: 



{n, j}(t, q) := J {1, v(p)} f(t, q, p) dp. 



It is well known that Eqs. (1.5, 1.3) amount to the existence of vector and scalar 
potentials such that 

B = curlA, E = —d t A - grad <I>. (1.6) 

The assumptions below the ID model are the following: all variables depend on 
only one space variable, denoted x, and the electrons are monokinetic in the direc- 
tions transversal to x. This is physically justifed by the fact that all the phenomena, 
especially heating, are much more rapid along the direction of propagation of the 
laser wave than in the transversal directions. So, the distribution function becomes: 

f(*,q,p) = f(t,x,p x )S(p± - p {t,x)) . (1.7) 

The function po(t,x) can be determined by Hamiltonian- mechanical considera- 
tions 16 . The Hamiltonian for one particle is H :— 7 m c 2 — e $, and (1.1) reads: 

g + [H,f]=0, 
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where [•,•] is the Poisson bracket. Then, by Hamilton's equation, the transversal 
component of the canonical conjugate momentum P c := p — e A is conserved: 

dP c = _c)H_ dP c± _ dH_ = 
dt <9q dt dq± 

By a suitable change of referential, we can suppose that P c j^ = 0; and by imposing 
the Coulomb gauge divA = 0, that A x = 0. Hence, po(t, x) — eA(t,x). 

For the sake of simplicity, we shall assume in this work that the pump wave is 
linearly polarised in a direction which we call y; however, the forthcoming compu- 
tations can be easily generalised to an arbitrary polarisation. Under these circum- 
stances, Eq. (1.6) becomes: 

B Z = 8 X Ay, Ey — —OtAy, E X = —d X $>, E Z — By = 0, 

which allows to recast the Vlasov equation (1.1) and the two remaining Maxwell 
equations (1.2, 1.4) as the following system: 

v + ^«_ t ( E , + i^)|U , (1 . 8) 

n 7 Ay (1.9) 



dt 7717 dx \ 7717 dx J dp. 



d Ay 2 d Ay 



dt 2 dx 2 77ie 

dE x e 



dt eo 



jx (1.10) 



dE x e 

— — = — (next - n) . (1.11) 
dx eo 

The Lorentz factor 7 and the density, quasi- density and flux are now given by: 



772 e 2 A 2 

7 = \/i + ^ + ^ri; (1-12) 

771 z C 771 z c 

{n,n 7 ,j x } = [ jl,-, — 1 fdp x . (1.13) 
J { 7 7717 J 

The equations (1.10) and (1.11), which are relative to the same variable E x , arc 
redundant under regularity conditions. Eq. (1.10) is a (simple) evolution equation, 
while (1.11) is interpreted as a constraint. As usual, the satisfaction of this constraint 
at t = implies its satisfaction at any time, thanks to (1.10) and to the continuity 
equation d t n + d x j x = 0. 

For a general polarisation, (1.9) would be duplicated, with a similar equation 
for A z . Quadratic terms in A y and A z would be added in (1.12) as well as in the 
third term in (1.8). The reader will convince himself that the study of this slightly 
more complicated system is completely similar to that of (1.8-1.13). From now on, 
we shall omit the subscripts in p x , E x , j x , A y . 

Let us notice that the model (1.8-1.13), as well as its extended version for a 
general polarisation, are classes of exact solutions to the relativistic Vlasov-Maxwell 
model, without any approximation. They are, as far as we know, the simplest exact 
solutions beyond ID Vlasov-Poisson models. 
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4 J. A. Carrillo and S. Labrunie 

1.1. Discussion of the relativistic character 

The model (1.8-1.13) features a strongly non-linear coupling between the kinetic 
and electromagnetic variables, through the Lorentz factor (1.12). This phenomenon 
makes this system difficult to study on the theoretical level, but also to solve nu- 
merically: no splitting between the variables is possible. 16 This is why two reduced 
models have been defined by physicists: 

(1) The non-relativistic model (hereafter denoted NR) approximates the relativistic 
dynamic by the Newtonian one by setting 7=1 everywhere. It is physically jus- 
tified when the temperature is low enough, so that the proportion of relativistic 
electrons is negligible, and the intensity of the pump wave is small. 



(2) The quasi-relativistic a model (QR) consists in approximating 7 by J 1 + 



in the second term in (1.8) and in the definition of j, and setting 7 = 1 in the 
third term in (1.8) and in the definition of n 7 , which amounts to setting n 7 = n. 
It is acceptable when the proportion of ultra-relativistic (v ~ c) electrons is 
negligible and the pump intensity is moderate. 

The original model, with 7 defined by (1.12) will be referred to as fully relativis- 
tic (FR). We remark that the NR model is a class of exact solutions to the non- 
relativistic Vlasov-Maxwell system, i.e. (1.1-1.5) with v(p) = p/m. By contrast, 
the QR model is only an approximation to the FR one. What makes its interest for 
the applied mathematician — besides its widespread use for simulation — is that it 
already contains certain features of higher-dimensional relativistic Vlasov-Maxwell 
systems, while being simpler to study. 

1.2. Rescaled equations and Cauchy problem 

The set of equations (1.8-1.13) can be simplified by introducing some rescaled 
variables. Let n be the unit of density; we choose the units for the independent 
variables as: 



and for the dependent variables as: 

— mc 2 n — mc — n 

E =\ > A = . /= > J = nc. 

V So e mc 

Keeping the same notations for the rescaled variables, we obtain the rescaled system: 





p = mc; 




(1.14) 



dE 
~dt 



J 



(1.15) 



a Also called semi-relativistic by some authors. 
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d 2 A d 2 A 

where the flux j and the quasi-density n 7 are defined as 

j(t,x):= I — f(t,x,p)dp, n 7 (t,x):= [ — f(t, x,p) dp. 
J 7i J 72 

Of course, suitable initial conditions are supplied, namely 

f(0,x,p) = f (x,p), A(0,x)=A (x), d t A(Q,x) = A {x). (1.17) 
For coherence, the initial electrostatic field must be given by: 

E (x) = E (0) - J (n (y) - n cxt (y)) dy, where: n (x) := J fo(x,p) dp. (1.18) 

This guarantees that the Poisson (or Gauss) equation 

dE 
dx 

will hold at any time provided the continuity equation is satisfied. 

As far as the relativistic character of the particles is concerned, the three versions 
of the model are respectively: 



= n ext -n, where: n(t,x) := J f(t,x,p)dp, (1-19) 



• NR: 7i = 72 = 1. 



QR: 7i = V 1 +P 2 , 72 = 1- 



• FR: 71 = 72 = y/l+jP + A 2 . 

Let us note that we shall investigate the existence of two classes of solutions 
for the system (1.14-1.16): periodic solutions, corresponding to initial data that are 
periodic in space, with a given period L, and "open-space" solutions, i.e. solutions 
of finite mass and energy. In both cases, we always assume that n cxt is at least 
bounded. Moreover, in the periodic setting, we assume that it is periodic and 

r-L 

(next (a) - n Q (x)) dx = 0, 



/ 

Jo 



10 

so that Eo is indeed periodic. 

To the best of our knowledge, this is the first mathematical work on this partic- 
ular Vlasov-Maxwell system (1.14-1.16): most of the previous mathematical works 
on reduced Vlasov-Maxwell models deal with systems living in two or one- and- a- half 
dimensions. 15,14,13,2 In our case, the lower dimension is, to some extent, compen- 
sated by a stronger nonlinearity. This system shares some common features with the 
Nordstrom- Vlasov system recently studied in 7,8 . Both systems present a Vlasov 
equation coupled to a wave equation whose right-hand side depends on the charge 
density. The main differences between them are the gravitational character of the 
Nordstrom- Vlasov system, our coupling with the Poisson equation and our more 
complicated right-hand side in the wave equation. 
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6 J. A. Carrillo and S. Labrunie 

In this article, we shall only investigate the existence and uniqueness of solutions 
to the NR and QR models. The differences between them can be bridged easily by 
using the shorthand notation 

p := p (NR case), p := — ^=^= (QR case). 

On the other hand, the method presented in this article cannot apply directly to 
the FR model, because of its much stronger and more non-linear coupling between 
the kinetic and electromagnetic variables. 

The next Section is devoted to reviewing some basic estimates for the Vlasov 
equation, and the different notions of solutions we will deal with. In Section 3 wc 
use the procedure of 9 to prove the global existence of weak solutions in the NR 
and QR cases, which are characterised as unique fixed points of a certain operator. 

This result is improved in two ways in Sections 4 and 5 respectively. We first 
prove the local- in-time existence of characteristic solutions in the NR case. The main 
difficulty in proving global existence of characteristic solutions consists in controlling 
the second space derivative of the vector potential A or, equivalently, the first space 
derivative of the density n. These difficulties are analogous to those found in the 
two and one-half dimensional Vlasov-Maxwell systems studied in 15 > u ,i3,2_ I n the 
QR case, we obtain the existence result thanks to a good integral representation of 
the second derivative of A reminiscent of similar ideas in 14 , used recently in 8 . 

Finally, the last section is devoted to obtaining a global energy functional for 
the three cases. It is conserved in time by characteristic solutions in the NR and 
QR cases. This fact, together with now standard relative entropy arguments, leads 
to the L p -nonlinear stability of a family of steady states in the periodic setting. 

2. Solutions to the forced Vlasov equation 

In this Section, we assume that the fields (E, A) are given and we introduce several 
notions of solution to (1.14) and summarise their associated regularity properties. 

Whatever their regularity, we are mainly interested in global solutions, i.e. which 
exist for any time. However, our estimates will generally not be uniform in time and 
thus, wc fix an arbitrary target time T and we look for solutions defined on [0,T], 
for all T > 0. We use the following notations for functional spaces: 

• C m (0,T; X): the space of m times continuously differentiable functions from 
[0,T] to the Banach space X; 

• VF fe >°°(R): the space of functions from R to R with all derivatives (in the sense 
of distributions) bounded, up to order k; 

• C b fe (R) = C k (R) n W k '°°{M.): the space of k times continuously differentiable 
functions with all derivatives bounded; 

• W£ ,tx> (R), C£(R): the subsp aces of functions which have the space period L. 

We will denote by ||-|| t the norm in L°°((0, t) x R), for t < T. 
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2.1. The characteristic system 

From now on, we denote 

F(t, x) := E(t, x) + A(t, x) d x A(t, x) (2.1) 

the force generated by the fields (E, A). The characteristic system associated to the 
transport equation (1.14) reads: 

dX ft , , dP 

(2.2) 



ds 

X(t) = x, 



= -F(s,X(s)), 

ds 

P(t) = p. 



Global existence and uniqueness of solution to the above system is ensured by 
assuming that the force field is continuous in time and globally Lipschitz in space; 
in turn, a sufficient condition for this is: 



AeC° (Q,T-W 2 '°° 



(2.3) 



Under this assumption, the unique solution to the characteristic system (2.2) de- 
noted by (X(s;t, x,p),P(s;t, x,p)) becomes (at least) a continuous function in all 
its variables. We shall also consider a stronger regularity condition 

(2.4) 



E e C° (0, T; C£(R)) , A E C° (0, T; C 2 ( 
From the uniqueness of the solution, we deduce a periodicity result: 

Lemma 2.1. If the force field is periodic in space, i.e.:V(t,x), F(t,x+L) = F(t,x), 
then the following identity holds for all s, t, x, p: 



X(s; t,x + L,p) = X(s; t, x,p) + L, P(s; t,x + L,p) = P(s; t, x,p). 



(2.5) 



The divergence of characteristics generated by different force fields is measured 
in a classical way (see for instance 9 or 2 and references therein) . 

Lemma 2.2. Let (X 1 ^ 1 ) and (X 2 ,P 2 ) be the characteristics associated to the 
respective forces F\, F 2 G C° (0, T; W 1 '°°{R)) . Then, the following inequalities hold 
for all (t, x,p): 

\X 1 {0;t,x,p)-X 2 {0:t,x,p)\<t f ||iq - F 2 || a ds, (2.6) 

Jo 

\P\0;t,x,p)-P 2 (0:t,x,p)\< f ||iq - F 2 \\ s ds. (2.7) 

Jo 

Proof. Eq. (2.7) clearly stems from the integration of the equation for P in (2.2). 
In the NR case, Eq. (2.6) immediately follows. In the QR case, we notice that 



_d 

dp 



l 



(1+ P 2) 



2^/2 



< 1, 



so 



p 2 



< IP 1 — P 2 \, hence (2.6) by integration. 



□ 
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In a very similar fashion, the derivatives of the characteristics can be estimated 
in terms of derivatives of the force field. 

Lemma 2.3. The solution (X,P) to (2.2), with F e C° (0, T; W 1 '°°(R)), is Lips- 
chitz w.r.t. the variables (x,p), and thus a.e. differentiable. The derivatives satisfy 
the following bound for a.e. (x,p) and all < r < t: 



max 



dX 



dX 

-dp- {T) 



dp 



op 

-dp- {T) 



< e 



(2.8) 



Proof. Consider two final conditions (xi,pi) and (x 2 ,p 2 ). We use the shorthand: 

(X < (*),P i (s)):=(X(s;i,a;i,pi),P(s;t,a;i ) p i )) ) i = 1, 2; 

and we denote by AF(s) the Lipschitz constant of F(s, •), i.e. sup x \d x F(s, x)\. 
By integrating the characteristic system (2.2), we get: 

\X 1 (t)-X 1 (t)\ + \P 1 (t)-P 2 (t)\ 

< \X1 - X 2 \ + \P1 ~P2\+ J {\Pl(s) - P 2 (S)\ + \F(8,X!(8)) - F(s,X 2 ( S ))\} ds 

< \xi-x 2 \ + \pi- P 2\+ J {\Pi(s)-P 2 (s)\+AF(s) \Xtis) - X 2 (s)\} ds 

< \xi - x 2 \ + \ Pl -p 2 \ + (1 + \\d x F\\ t ) JjlX^s) - X 2 (s)\ + \Piis) - P 2 (s)\} ds. 

Hence, by Gronwall's lemma: 

|Xi(t) - X 1 (t)\ + |P!(r) - P 2 (t)\ <(\ Xl - x 2 \ + \ Pl - p 2 \) e<*" T > (i+ll^^llO, 

which proves the Lipschitz character of the functions X, P and the quantitative 
estimate (2.8). □ 

2.2. Characteristic, mild and weak solutions 

In this work, we shall always assume the following two hypotheses about the initial 
distribution function /o. 

Hypothesis 1. In the open-space case, fo € L 1 (K 2 ); in the periodic case, 

f G Ll(R 2 ) := {/ e Ll c (R 2 ) : f(x + L,p) = f(x,p) a.e. m (x,p) G M 2 

^d /| (Oii)xR GL 1 ((0,L)xR)|. 
In both cases, (d x f ,d p f ) G Lj^R 2 ). 

Hypothesis 2. TTiere exists a continuous, positive, even function g(p), which more- 
over is decreasing in \p\ and satisfies 

\p\ g(p) dp < oo, 
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such that 

fo(x,p) < g(p), \d x f {x,p)\ < g(p), \d p f Q (x,p)\ < g(p). 

Let us notice that these hypotheses imply, in the open-space case, that /o £ 
W^{R 2 ) n M /1 ' 00 (R 2 ); the initial density n £ L 1 (R) n M /1 ' <X> (R); and the initial 
electrostatic field, given by (1.11), E £ W 1 ' 00 ^) if n cxt £ £ £ C£(R) if 

«cxt £ C°(R). In the periodic case, one has similar properties in terms of periodic 
spaces; indeed, electrical neutrality ensures that n—n cxt admits a periodic primitive. 

For further reference, let us note that the function g r (p) defined as 

9r(p) = .9(0) for \p\ < r, g r (p) = g(\p\ - r) for \p\ > r, 

satisfies 

f g T {p) dp = 2rg(0) + f g(p) dp. (2.9) 
Jm Jr. 

The first notion of solutions which we consider is the one given by the usual 

characteristic method. 

Definition 2.1. Given a force field F £ C° (0, T; W 1>00 (R)) and / satisfying Hy- 
potheses 1-2, we define the characteristic solution of the Vlasov equation (1.14) 
as: 

f(t,x,p) := f o {X(0;t,x,p),P(0;t,x,p)), (2.10) 

where (X(s; t, x,p), P(s; t, x,p)) is the unique solution to (2.2). Moreover, assuming 
F £ C° (0, T; Cl(R)) and f £ C^R 2 ), then we shall often refer to the characteristic 
solution as a classical solution. 

The previous definition gets clarified by the following result. 

Lemma 2.4. Under the assumption (2.3), the characteristic solution f belongs 
to W 1 ' 00 (0, T; W 1 '°°(M 2 )) . If, moreover, f £ C\R 2 ) and (24) holds, then f £ 
C 1 (0,T;C^(M 2 )) and (1.14) « s satisfied in the classical sense. 

Proof. According to classical dynamical system theory, (X(s; t, x,p), P(s; t, x,p)) 
is C 1 w.r.t. s, and also w.r.t. t given the symmetry of these variables. The Lipschitz 
character in (x,p) was obtained in Lemma 2.3. The first part of the conclusion then 
follows from Hypothesis 2, and the fact that the composition of Lipschitz functions 
is Lipschitz. 

Moreover, under (2.4), the force is C 1 ; as p *—> p is C°°, we deduce that the 
solution to (2.2) is C 1 . The last statement then follows from the chain rule. □ 

Mild solutions are introduced for relaxing the assumption of differentiability of 
the characteristics but keeping the fact that they define a family of Lipschitz home- 
omorphisms in phase space (see 1 and references therein). In fact, this definition 
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can be rephrased using the concept of push-forward of a density through a map, 
which is quite well-known in mass transport theory. 18 The pushed- forward measure 
T#p of a given measure p in R" assigns mass 

T #P [K] := p[T-\K)} 

to each Borel set K C R™. By this property, it satisfies that 

/ V>d(T #/9 )= f (^oT)dp 
for all test functions tp £ C°(R n ). 

Definition 2.2. Given a force field F G C° (0, T; VF^°°(R)) and / G ^(R 2 ), 
we say that a weakly continuous function f(t,x,p) G C w ([0, T]; L\ oc {U?)) is a miZd 
solution of the Vlasov equation (1.14) if it satisfies 

/ ip(x,p)f(t,x,p)dxdp = / ip{X(t;0,x,p),P{t;0,x,p)) f (x,p)dxdp (2.11) 
Jr 2 Jm. 2 

for all test functions ip G C°(M 2 ) and all t > 0, i.e., 

f(t,x,p) = {X(t;0,x,p),P{t;0,x,p))#f 

for all t > 0. 

Taking into account the change of variables formula for Lipschitz functions, 11 wc 
deduce that a characteristic solution is a mild solution of the Vlasov equation (1.14). 

We can relax even more the assumptions on the force field and talk about dis- 
tributional solutions for the Vlasov equation (1.14). 

Definition 2.3. Given a force field F G L°°((0,T) x R) and / G i^QR 2 ), wc 
say that f(t,x,p) G Ll oc ((0,T) x R 2 ) is a distributional solution of the Vlasov 
equation (1.14) if it satisfies 

-[SM +id i- F %) Iix dpdt= L m x - p), ° dxdp <2 - 12) 

for all test functions ip G C c °°([0,T) x R 2 ). 

It is easy to check that any mild solution is a distributional solution. Moreover, 
in the particular case of characteristic solutions, in which / G T / F 1, °°(R 2 ), / belongs 
to W 1 ' 00 (0,T; W 1 -°°(R 2 )) by Lemma 2.4. Therefore, one can check from the weak 
formulation (2.12) that / satisfies the Vlasov equation (1.14) as an equality almost 
everywhere of locally bounded functions on (0, oo) x R 2 . 

2.3. A priori estimates 

Since bounds for functions in C^J(R) or W /1 '°°(R) can often be obtained in the same 
manner, we shall treat in the following the two types of characteristic solutions 
(classical and mild) together. 
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Thanks again to Hypothesis 2, one can define the density and flux: 

{n,j}(t,x) := f {l,p}f(t,x,p)dp. (2.13) 

JR 

Indeed, we have the following more general lemma for moments of characteristic 
solutions of (1.14). Let us denote by rrik(t,x) the moment of order k 6 N of a 
solution f(t, x,p) with a given force field F. 

Lemma 2.5. Let f be a characteristic solution of the Vlasov equation (1.14) w ^ 
force field F E C° (0, T; W 1,00 (K)) . Assume that the function g(p) from Hypothesis 2 
has the moment of order k bounded, i.e., 

M k ■= / \p\ k g(p)d P < oo, 

JR 

then rrik(t,x) is well-defined for all (t,x) and 

\m k (t,x)\ <M k + R k (M + M k ,t\\F\\ t ) (2.14) 
where R k (a,b) will be defined below. In particular, for k = we obtain 

\\n\\ t <M + 2g(0)t\\F\\ t . (2.15) 

Proof. From (2.10) and Hypothesis 2, we deduce: 

\p\ k f(t,x,p) < \p\ k g(\P(0;t,x,p)\) . 

Now, Eq. (2.7) with one of the force fields replaced by yields |P(0; t, x,p) — p\ < 
t \\F\\ t ; as g is decreasing, this gives: g (\P(0;t,x,p)\) < gt\\F\\ t (p), an d thus 

\p\ k f{t lX ,p)<\p\ k g tmt {p) 

which is clearly integrable in p. This proves that m k (t,x) is well-defined and 



\m k (t,x)\ < J \p\ k g t \\F\\ t (p) d P- 



Moreover, we notice that for k > 1 we have 

k 



[ \p\ k 9r( P )dp = Y^r k+1 +J2 C l rk ~ i I \p\ l 9(p)d P 

] l r ^+(j2cir k A [(l + \p\ k )g(p)dp 
\i=i ) Jr 



" k + 
2 5(0) 



< 



k + 

:=R k {M + M k ,r). (2.16) 

Finally, combining previous inequalities, we obtain (2.14). The estimate on the 
density (2.15) follows directly from previous arguments and (2.9) which defines the 
function Rq. □ 
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We can also estimate the divergence of moments corresponding to two different 
solutions of the Vlasov equation (1.14). 

Lemma 2.6. Let /i, f 2 be characteristic solutions of the Vlasov equation (1.14) 
with forces F\, F 2 <G C° (0,T; W 1,00 (R)) respectively. Assume that the function g(p) 
from Hypothesis 2 has the moment of order k bounded, then 

\m 1 , k -m 2 ,k\{t,x) <R k (M + M k ,max.\\Fi\\ t } J - F 2 \\ s ds. (2.17) 



Proof. To estimate m\ t k — rn 2 ,k, we consider the characteristics (A^P 1 ) and 
(X 2 ,P 2 ) associated to F\ and F 2 respectively. Using the shorthand notation for 
the characteristics (X^Pj) = (X*(0; t, x,p), P l (0; t, x,p)) , we write: 



(mi, fc - m 2 , fc ) (t, a:) = / |p| fc [/o (Xq 1 . Pq 1 ) - / (X 2 , P 2 )] dp 



and thus, 



|mi, fc - m 2 ,fe| (i,a;) < / |p| fc 



EM) 



+ 



<9p 



X 2 P 



l-^o — I 



I Po-Po I 



dp, 



where we have made use twice of the one-dimensional Taylor-Lagrange formula; X, 
respectively P, lie between X^ and X 2 , resp. Pq 1 and Pq. Then we invoke (2.6-2.7) 
to bound 



\mi t k - m 2 ,fe| (M) < 



jf|p| fc {s(|^o|)+*s(|£|)} dpx| 11^1-^11,(1*. (2.18) 

Applying again Eq. (2.7) with one of the force fields replaced by yields | Pq — p | < 
t ||P;|| t ; as g is decreasing, this gives: g (|Po|) < gt\\Fi\\ t {p), an d then: 

ff(|P|) < maxj^dPo 1 !) ,.g(|P 2 |)} < maxg t m \ lt (p). 

Moreover, r \— ► g r (p) is an increasing function of r, for all p; this implies: 

/ \p\ k {9(\Po\)+tg(\P\) \ dp< (1+t) / |p| fe max fft || Fi || (p) dp 

JR L JR 1 

< (1+t) max / |p| fc 3t || Fi || t (p)dp. 

Inequality (2.17) is obtained by using the last lines in Lemma 2.5. □ 

A small variation on the above arguments allows us to prove that in fact, the 
density and flux are regular enough to satisfy the continuity equation. This is impor- 
tant in order to be able to say that the description via the Poisson equation (1.11) 
is equivalent to the Ampere equation (1.10). 
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Corollary 2.1. Under Hypotheses 1-2 and the assumption (2.3), resp. (2.4), n and 
j are Lipschitz, resp. continuously differentiable in space and time and they satisfy 



dn dj 
h — = 0. 

dt dx 



(2.19) 



Moreover, there is a constant C, depending only on the majorising function g(p), 
such that: 

dn 



dx 



< Ce'( 1+l|S:cF|lt ). 



(2.20) 



Proof. When the solution is classical, we can apply the chain rule and estimate 
directly the derivative in x of the solution f(t, x,p) to get 



df 



dx 



< 



sup 



dfo 



dx 



( sup 



dP 



dfo 



dpo 



Lemma 2.3 and Hypothesis 2 imply that (see previous Lemma): 



df 



dx 



<eKW^ w ,( p ), 



and thus d x f is integrable. Moreover, Lebesgue's dominated convergence theorem 
implies that n is Lipschitz or differentiable with respect to x and hence (2.20). 
Similar arguments apply to dtn and d x j and the continuity equation (2.19) becomes 
an easy consequence of Eq. (1.14) upon integration on p and using Hypotheses 1-2 
on the initial data. In case the solution is only mild under assumption (2.3), one can 
reproduce the Lipschitz bounds by estimating the difference n(t, Xi) — n(t, x 2 ) using 
analogous arguments to previous Lemma 2.6, we leave the details to the reader. □ 

And it stems from Lemma 2.1 that: 
Lemma 2.7. Under the hypotheses: 

V(x,p), f {x + L,p) = f (x,p); V(t,x), {E,A}(t,x + L) = {E,A}(t,x); 
there holds: 

V(i,x,p), f(t,x + L,p) = f(t,x,p), {n,j}(t,x + L) = {n,j}(t,x). 



Let us finally remark that thanks to the mass conservation property of the 
Vlasov equation, we estimate the integrals of n. In the open-space case, we have 

rb 

n(x)dx<M, (2.21) 



/ 

J a 



where M is the total mass of /q. In the periodic case, 



n(x) dx < M 



(2.22) 
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where M is now the mass of /o over one period, and \r~\ is the smallest integer 
larger or equal to r. 

Moreover, in the open-space case we can estimate moments in x under suitable 
additional assumptions on the initial data. In fact, let us consider the following 

Hypothesis 3. fo(x,p) < g(x) g(p) in R 2 . 

By following the same lines of argument as in Lemma 2.5, we can prove: 

Lemma 2.8. Let f be a characteristic solution of the Vlasov equation (1-14) in the 
open-space case with force field F <E C° (0, T; M /1 '°°(R)) and assume that Hypoth- 
esis 3 is satisfied. Then, there exists a constant Z , depending polynomially on the 
first two moments of g(p), t and \\F\\ t , such that 



This can be generalised to all moments of / and its derivatives, under 

Hypothesis 4. \d x fa(x,p)\ < g(x) g(p) and \d p f (x,p)\ < g(x) g(p) in R 2 . 

Lemma 2.9. Assume that Hypotheses 3 and 4 hold, and that g has its moment of 
order m bounded. Then, for all k, £ < m, there exists a constant Z^.t such that, for 
all t G [0,T]: 



and a similar bound holds for d p f. 

3. Iterative procedure and global weak solutions 

In this section, we present an iterative procedure to solve the ID Vlasov-Maxwell 
system (1.14-1.16) for the NR and QR cases. 

First, we define the iterative procedure based on the Cooper-Klimas 9 approach. 
Then, we will derive estimates on the fields that allow us to obtain a limit by tele- 
scopic series. However, these estimates will not allow us to get a global characteristic 
solution, due to the lack of a global-in-time estimate on the space derivative of the 
density, which is needed to control that of the right-hand side in (1.16), and thus the 
second space derivative of the vector potential. Therefore, at this level of generality 
we are only able to obtain global weak solutions. Improvements of this basic result, 
namely local (in the NR case) and global (in the QR case) existence of classical and 
mild solutions will be postponed to the next two sections. 




(2.23) 
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We now fix an initial condition / for the distribution function, as well as initial 
data (A ,A ) for the vector potential, no and E are the density and electrostatic 
field given by /o, as defined in (1.18). Let us remind that we always assume that 
the external density verifies n cxt G L°°(R). We also fix a target time T, and, if wc 
are interested in periodic solutions, a space period L. 

3.1. Definition of the recurrence operator 

Given the field pair (E,A) G C° (0, T; W 1 '°°(R) x W 2 '°°(R)) , one constructs 
(E',A') =C(E,A) as follows. 

(1) The characteristic system (2.2) is solved, with the force F given by (2.1). 

(2) One computes the characteristic solution to Vlasov's equation by (2.10) and its 
density and flux (n,j) by (2.13). 

(3) Finally, E' and A 1 arc computed as: 

E'(t,x):=E {x)+ [ j(s,x)ds, (3.1) 
Jo 

1 f f x+t ■ 

A'(t, x) := - U (x + t) + A (x -t) + J A (y) dy 

ft px-\-t — S -\ 

^ (nA)(s,y)dyds\. (3.2) 



IT 

JO Jx+s 



From Corollary 2.1 and the Duhamel formulae (A.2-A.4), we immediately deduce: 

Theorem 3.1. If (A Q ,A ) G VF 2 '°°(R) x W 1 '°°(M.), the operator C maps 
C° (0, T; W 1 '°°(R) x W 2 -°°{R)) to itself. If(A ,A ) G C b 2 (M) x cf(R), f G C X (R 2 ), 
and n cxt G C${R), then C maps C° (0,T;C^(R) x C b 2 (M)) to itee//. 

Moreover, the Poisson equation (1.19) is satisfied for the pair (E',n); and the 
time derivative of A' is bounded m x, i.e. A' G C 1 (0, T; L°°(R)) or C 1 (0, T; C{?(R)) , 
euen i/ A rfoes not belong a priori to smc/i a space. 

This ensures that the operator £ can be iterated. Of course, we shall need some 
quantitative estimates. To establish them is the goal of the next Subsection. 
Moreover, we deduce from Lemma 2.7 the following 

Corollary 3.1. If (A ,A ) G W^'°°(R) x W^'°°(R), the operator C m,aps 
C° (o,T;Wl'°°(R) x W| ,0C (R)) to iteetf. //(A ,i ) G C£(R)xCi(R), fo G Ci(R 2 ) 7 
and n oxt G C° L {R), then C maps C° (0,T;Cl(R) x C£(R)) to itself. 



3.2. A priori estimates 

In the sequel, we shall always assume at least that (A , A ) G VK 2 '°°(M) x W 1,00 (M). 
The constants denoted C or Cj may vary from one line to the next, and depend on 
the initial conditions, T and L (but on nothing else). 



February 2, 2008 8:39 WSPC/INSTRUCTION FILE carlabnew 



16 J. A. Carrillo and S. Labrunie 



First, we fix (E,A) G C° (0, T; W 1,0 °(R) x W 2 '°° 
F is the force field corresponding to (E, A). 

For the electrostatic field we have the following properties: 

Lemma 3.1. The following estimates hold: 

\\E'\\ t <C + Ci fnds; 
Jo 

<Co + Ci||F|| t . 



and {E',A') := C(E,A). 



BE' 



dx 



(3.3) 
(3.4) 



Proof. To obtain the first estimate, we proceed as in 1 by duality. Given ip e L 1 (l 
we apply changes of variables to get 

/ / j(s,x)<p(x)dsdx = / f (x,p) / P(s;0,x,p)<p(X(s;0,x,p))dsdxdp 
JrJo Jr 2 Jo 



rX(t;0,x,p) 

fo{x,p) I ip(u)dudxdp, 



and thus, using Hypothesis 2 



j(s, x)ip(x) ds dx 
r Jo 



rX{t-fl,x,p) 



ip(u) du 



da; dp 



< 



/ fo(x,p) / 
Jr 2 Jx 

/ g(p) / \f( u )\ / x(u;S(t,x,p))dxdudp 
Jr Jr Jr 



where x( u : S(t, x,p)) is the characteristic function of the set S(t, x,p) = {\u — x\ < 
\X(t; 0, x, p) — x\ } as a function of u. Using the bound on the divergence of forward- 
in-time characteristics as in Lemma 2.6, we obtain 

\X(t;0,x,p) - x\ < t [ \\F\\ s ds and f x(u; S(t, x,p)) dx < 2t [ ||F|| s ds; 

Jo Jr Jo 

and therefore, we deduce 



Jo 



j(s,x) ds 



< It 



[ g(p)dp [ 

Jr Jo 



\F\\ s ds. 



Finally, (3.3) is deduced directly from (3.1). 

Now, Corollary 2.1 assures that the continuity equation (2.19) is satisfied. Taking 
into account that equation and the Ampere equation (3.1), then 

dE' 



dx 



= n ex t - n, 



holds. Estimate (3.4) follows from (2.15) and n oxt £ L c 



□ 



For the vector potential, we have a similar result. 
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Lemma 3.2. The following estimates hold: 

(3.5) 

Jo 

dA' 



A'|| t <C + Ci / \\A\\ s ds, 
Jo 

<Co + \\A\\ t (Ci+Ci J ||F|| s ds) , 
<C + \\A\\ t (Ci + C* J*\\F\\ s ds^ . 



dx 
dA' 



(3.6) 
(3.7) 



Proof. From (3.2), we bound 

\A'(t,x)\<C + /Vll- ( f X+t S n(s iy )dy] ds. 

J0 \Jx+s-t / 

Using (2.21) or (2.22), we then bound the integral of n by M or M\2T/L]. 
Hence (3.5). 

Let us estimate the derivatives of A. d x A' is given by the Duhamcl formula (A. 3), 
with / replaced by — n A. Hence the majoration: 



dA' . 



<C + ||A|| t / (n(s,x + s-t)+n(s,x-s + t))ds. 



Now, using the uniform bound on the density (2.15), we obtain 

( (n(s,x + s-t)+n(s,x-s + t))ds<Ci+C 2 [ \\F\\ s ds. 
Jo Jo 

This gives (3.6). As the Duhamcl formula (A. 2) for the time derivative is very similar 

to (A. 3), one can establish (3.7) by the same reasoning. □ 

Now, we consider two field pairs (£i, A\) and (^2,^2), with the corresponding 
forces Fi and F 2 , and we set (E[, A[) := C{E U A{), (E' 2 ,A' 2 ) := C(E 2 ,A 2 ). 

Lemma 3.3. The following estimates hold for all field pairs (Ei, A\) and (E 2 ,A 2 ): 
||^-^|| t <i2i(Mo + Mi,max||Fi|| t ) J ^ - F 2 \\ s ds, (3.8) 

\\A[ -A' 2 \\ t < \\A 2 \\ t (C + d max ||F|| t ) jT* - F 2 \\ s ds 

L 



+ c 2 



\A2-AAds. 



(3.9) 



Proof. From (3.1) we deduce 

\\E[-E 2 \\ t < f \j 1 (s,x)-j 2 (s,x)\ds, 
Jo 

and thus (3.8) is a simple consequence of Lemma 2.6. 
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From (3.2), we have 

2(A' 1 -A' 2 )(t,x) = / (n 2 A 2 -n 1 A 1 )(s,y)dyds. 

JO Jx+s-t 

Writing: n 2 A 2 — n\ A\ = (n 2 — n\) A 2 — (A 2 — Ax) m, we bound: 

2 Pi -411* < / P2|| 8 / \n 2 (s,y)-n 1 (s,y)\ dyds 

JO Jx+s-t 

+ / — -4.1 || s / n 1 (s,y)dyds. 

JO Jx+s-t 

Using (2.21) or (2.22), the second line ol this inequality is easily bounded as 
AM f* \\A 2 -A 1 \\ s ds, with A = 1 or \2T/L]. Then, the first line is bounded 
thanks to Lemma 2.6: 

\ ni -n 2 \(t,x) < (Co + d max||F4) J ||iq - F 2 || s ds, (3.10) 



and finally: 



K - A' 2 \\ t < \\A 2 \\ t (C + Ci max ||F,|| t ) J* ||^ - F 2 \\ s ds 



+ 



AM 



Jo 



2 - Ai|| ds, 



(3.11) 



which is (3.9). 

Now, we estimate the difference of the derivatives. 
Lemma 3.4. There holds, for any (E\,A\) and (E 2 ,A 2 ): 

< (Co + Cj max||F,|| 4 ) J \\F 2 - F^ ds, 

< (Co + Ci max||^|| t ) ^ \\A 2 - A^ ds 
+ \\A 2 \\ t J*\\F 1 -F 2 \\ a ds}, (3.12) 

< (Co + d max||fi|| t ) jjf \\A 2 - ds 

+ ||A 2 || t |j| J F 1 - J F 2 || s ds|. (3.13) 

Proof. The Poisson equation and (3.10) imply (3.11). From the Duhamel for- 
mula (A. 3), with / replaced successively with — n\ A\ and — n 2 A 2 , we derive: 



dE[ 


8E' 2 


dx 


dx 


dA[ 


dA' 2 


dx 


dx 


dA[ 


dA' 2 


dt 


m 



dx dx 



dA[ dA' 2 \ 

z \(t,x) = J (n 2 A 2 -mAi)(s,x + s-t)ds 

- / (n 2 A 2 - m Ai)(s, x + t - s) ds ; 
Jo 
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once more, we write n 2 A 2 — n\ A\ = (n 2 — ni) A 2 — (A 2 — A\) m, which gives 

\\A4 t fwrn-nzlds+WmWt f \\A, - A 2 \\ s ds. 
Jo Jo 

Using the bounds (3.10) for the first term, and (2.15) for the second yields 
-q£--q£ <\\M t (Co + C.m^mi) jf -F 2 \\ s ds 

+ (C 2 ||*i ||t + C 3 ) /"* ||^i - A 2 || s ds, 
Jo 

which implies (3.12). Once again, the similarity ol the formulae (A. 2) and (A. 3) 
allows to deduce (3.13). □ 



3.3. Convergence of successive approximations 

We start from the initial data (E (x), A (x)), by extending them to constant-in- 
time functions over [0, T] x R. Then, we construct a sequence (E k , A k ) k& ^ by the 
recurrence formula: 

(E k+1 ,A k+1 ):=£(E k ,A k ), Vfc > 0. 

Of course, we set F k := E k + A k d x A k . 

To establish the convergence of (E k ,A k ) keN , the following result will be useful. 
It is easily proved by induction. 

Lemma 3.5. Let (ttfc)fceN be a sequence of positive functions, u k : [0,T] — > R+, 
satisfying 

(i) Vte[0,T], u {t)<c, 

(ii) Mk G N, Vi G [0,T], u k+1 (t)<a + b u fe (s)ds, 

Jo 

for some constants a, b, c G R + . Then the following estimate holds: 
Vfc G N, Vt G [0, T], Wt) <aV -T-+c- 



i=l 



fc! 



Zience the whole sequence is uniformly bounded by a constant u* which depends on 
a, b, c and T . 

Theorem 3.2. The sequence (E k , A k ) ke f$ converges uniformly in t and x towards 
a limit (E, A). 

(1) If(A ,A ) G W 2 ' co {R) x W 1 '°°(R), the sequences (E k ) and (A k ) converge re- 
spectively within C° (0, T; W 1 '°°(R)) and C 1 (0, T; L°°(R))nC (0, T; W 1 >°°(R)) . 

(2) If{A ,A ) G C 2 (R) x Cl(R), f G C X (R 2 ), and n cxt G C%{R), then the conver- 
gences take place within C° (0, T; C^R)) and C 1 (0, T; C£(R)) nC° (0, T; C£( 
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(3) The above conclusions are valid when all spaces are replaced by their periodic 
counterparts. 



Proof. First, Eq. (3.5) gives: 

Pfe+i|| t <C + Ci / ||A fc || 8 ds 
Jo 

and thus, the previous lemma ensures that ||^4fe|| t < A* uniformly in k and t. Then, 
Eqs. (3.3) and (3.6) imply: 

\\F k+1 \\ t <\\E k+1 \\ t + A* \\d x A k+1 \\ t < C 2 + C 3 f \\F k \\ s ds, 

Jo 

so, once more, ||-Ffc|| t < F». Using again (3.3) and (3.6) shows that the sequences 
||i?fe|| t and ||9 x j4fe|| 4 are uniformly bounded in k and t. 
All these bounds allow to rewrite (3.9) and (3.12) as 

Pfe+i -A k \\ t < C\ f \\F k -F k _ 1 \\ s ds + C 2 [ \\A k -A k -!\\ a ds; (3.14) 
Jo Jo 

d -^r- d -^T <^ j\\F k -F k ^\\ s ds + C 2 j\\A k -A k ^\\ s ds. (3.15) 
Then, writing 

F k+ i - F k = E k+ i - E k + A k+ i (d x A k+ i - d x A k ) + d x A k (A k+1 - A k ) , 
one deduces from (3.8), (3.14) and (3.15) 

\\Fw-F k \\ t <Ci f \\Fk-Fk-4. ds + C 2 [ \\A k -. 

Jo Jo 



A k - 



ds. 



(3.16) 



Hence, the sequence u k (t) := \\A k — A k -i\\ t + \\F k — F k _\\\ t satisfies: u k+ i(t) < 
b f*u k (s)ds, for some constant b. Lemma 3.5 with a = then implies u k (t) < 
b k T k /k\, i.e.: 



\A k+1 - A k \\ t < — — , \\F k+ i - F k \\ t < 



~fcT : 



so by (3.8) and (3.15) 



\E, 



fc+i 



E k L<C 



~k~T' 



dA k+1 dA k 



dx 



dx 



< C 



jyk rpk 

k\ ' 



As a consequence, the four sequences {E k ) km , (A k ) keN , (d x A k ) keN , (F k ) keN , 



all 



converge uniformly on [0, T] x M. Let E, A, B, F, be the limits; clearly B = d x A 
in the sense of distributions, and F = E + A d x A. Under the hypothesis (A , A a ) E 
W 2 '°°(R) x W 1 ' 00 ^), each term of the four sequences is in C° (0, T; L°°(R)), and so 
arc the limits. So: E G C° (0, T; L°°(R)) and A G C° (0,T; W 1,0O (R)). If (A ,A ) G 
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and n oxt G C£(R), the terms and the limits arc in C° (0, T; C°(R)) , i.e. £ G 
C° (0, T; C b °(R)) and A G C° (0, T; C£(R)) . 

With all these results, Eq. (3.13) shows that: \\d t A k+1 - d t A k \\ < Cb k T k /k\; 
since Aq is independent of time, the sequence (dtAk) k ^ also converges uniformly in 
C° (0,T;L°°(R)), resp. C° (0, T; C°(K)) , towards a limit which is necessarily equal 
to d t A. Thus, A eC 1 (0,T;L°°(R)), resp. C 1 (0,T;C°(R)) . 

Similarly, Eq. (3.11) implies that (d x Ek) ke j i is a Cauchy sequence and thus 
converges uniformly toward a limit which is necessarily equal to d x E. Hence, E G 
C° (0, T; W^°°(R)) or C° (0, T; C*(R)) . 

The last point easily follows from Corollary 3.1. □ 

3.4. Global existence and properties of weak solutions 

If (E,A) G C° (0, T; VF 1,00 (R)) x C° (0, T; TT /2,00 (M)) , then it is a fixed point of 
C within this space. Thus, the triple (f,E,A), where the function / is defined 
by (2.10), is a characteristic solution to (1.14-1.19). 

This fixed point and its associated characteristic solution, if they exist, are 
unique. Indeed, let (f\,E\,A\) and (f 2 , E2, A2) be two such solutions, with respec- 
tive force fields F\ and F 2 . Reasoning as in the proof of the above theorem shows: 

\\A! - A 2 \\ t + \\F t - F 2 \\ t < C f {\\A X -A 2 \\ s + \\F X F 2 \\ s } ds ; 

Jo 

hence A\ = A 2 and F\ — F 2 by Gronwall's lemma, then E\ — E 2 by (3.8), and 
finally fi = f 2 by the Cauchy-Lipschitz theorem. 

Without further a priori estimates on the second derivatives of the vector po- 
tential A, we will not be able to obtain global existence of characteristic solu- 
tions. However, the bounds in Subsection 3.2 show that the norm of C{E, A) in 
C° (0,T;L°°(R) x W^°°(R)), resp. C° (0,T;C£(R) x C^(R)), is controlled by the 
norm of (E, A) in the same space. These considerations suggest to extend the oper- 
ator £ to a bigger space in order to obtain weak solutions, which would also enjoy 
a uniqueness property. This programme cannot be achieved within the framework 
of C° (0, T; i°°(R) x W hoo (R)). The reason is the well-known lack of density of 
smooth functions within L°°(R) or T4 /1,00 (IR). But this obstruction is removed if 
one works within C° (0,T;C£(R) x C£(R)). 

Lemma 3.6. Assume (A ,A ) G C^(R) x C£(R), f G C^R 2 ), and 
n e xt G C°(R). The operator C can be extended to a continous operator from 
C° (0,T;C£(R) x Cl(R)) to itself, which satisfies the estimates of Subsection 3.2. If 
the initial conditions are periodic, then C maps C° (0,T;C£(R) x C^(R)) to itself. 

Proof. C is defined on the dense subspace Y := C° (0,T;C^(R) x C h 2 (M)) of X := 
C° (0,T;C h °(R) x Cl(R)), with values in X. From (3.8, 3.9, 3.12), we see that it is 
uniformly continuous, in the norm of X, on any set K n Y, where K is a bounded 
set of X. Hence, it admits a unique continuous extension from X to itself. 
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The extension procedure preserves the estimates of Subsection 3.2, as their r.h.s. 
are uniformly continuous in the norm of X on any K n Y. Finally, the closedness 
of C° (0, T; C£(R) x C^(R)) within X guarantees the invariance of this space by the 
extended operator. □ 

Hence, the following theorem, whose proof rephrases that of Theorem 3.2. 

Theorem 3.3. Assume (A ,A ) E C 2 (R) x C£(R), / E C^R 2 ), anrfn cxt e C£(R). 
XTie operator C admits a unique fixed point (E, A) e C° (0, T; C£(R) x C^R)) , *cft 
moreover belongs to C° (0,T;C£(R)) x C 1 (0, T; C b °(R)) . J/tfte imfwi conditions are 
periodic, so is the fixed point. 

Let us now check that this unique fixed point defines a distribution function / in 
such a way that the triple (/, E, A) is a global weak solution of the Vlasov-Maxwell 
system (1.14-1.19). 

Theorem 3.4. Under the hypotheses of Theorem 3.3, let (Ek,Ak)k a sequence con- 
verging to the fixed point (E, A); without loss of generality, we can assume that its 
terms belong to C° (0, T; C^(R) x C 2 (R)) . Let (fk)k be the associated sequence of dis- 
tribution functions obtained by the method of characteristics. Then, (fk)k converges 
uniformly in all its variables toward a function f E C° (0, T; C° (R 2 )), which does not 
depend on the sequence (Ek,Ak)k, and the triple (f,E,A) satisfies (1.14)-(1.19) in 
the sense of distributions. 

Proof. Repeating the argument of Lemma 2.6, we easily obtain the estimate 

\fk(t,x,p)-f e (t,x,p)\<(l+t)ma X g tllFillt (p) f \\F k - F t \\ s da, (3.17) 

i— k,t Jq 

from which follows that (fk{t)) k is a Cauchy sequence in L°° norm, uniformly in t. 
A similar argument shows that, if we consider another approximating sequence 
denoted by tildes, fk(t) — fk(t) will converge to zero, uniformly in x, p and t. Let 
us then call / the common limit of all the sequences (fk)k'i continuity in all 
variables follows from that of the f k ■ 

As fk is a classical solution to the Vlasov equation, it is also a distributional 
solution, i.e. an equation similar to (2.12) holds, with / and F replaced resp. with 
fk and Fk- For any test function ip, the integrals in this formula are taken over 
a compact subset of [0, T) x M 2 . Thus, the uniform convergence in i, p and t of 
the sequences (/fc) fe and (Fk) k ensures that (2.12) will also hold at the limit: the 
limiting triple (f,E,A) satisfies (1.14-1.19) in the sense of distributions. □ 

The above result can be improved in the following two ways. First, we show that 
f(t) will retain the mass of the initial distribution f . 

Proposition 3.1. In the periodic case, fk(t) converges toward f(t) in L^(K 2 ), 
uniformly in t, and the mass of f(t) over one space period is equal to that of fo- 
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In the open-space case, under Hypothesis 3, f k (t) converges weakly toward /(f) 
in L 1 (M 2 ) ; and the mass of f{t) is equal to that of Jo- 
Proof. Inequality (3.17) shows that fk(t, x, •) converges to f(t, x, •) in uni- 
formly in x and t. Hence the convergence of f k (t) toward f{t) in L 1 ((a, b) x K) for 
any a < b. 

Yet, the fk are classical solutions to the Vlasov equation. In the periodic case, 
they conserve the mass of /o over one space period. So, /(f) retains this mass. 

In the open-space case, taking into account the L 1 convergence in bounded 
intervals in space, it suffices to show that mass does not escape at infinity, uniformly 
in t. Using Lemma 2.8, we obtain: 

/ \f k (t,x,p)\dxdp < 4 / (\x\ + \p\)fk(t,x,p)dxdp<^->0, 

J\x\ + \p\>R H J\x\ + \p\>R & 

uniformly in k and f when R — > oo. Thus, it is a standard argument to check that 
J R 2 f{t, x,p)dxdp = lim fe J* R2 f k (f , x, p) dx dp = M . □ 

The second useful precision is that /(f) is "almost a mild solution" to (1.14). 

Theorem 3.5. The characteristics given by the Fk converge uniformly in all their 
variables toward a limit (X, P), and Eq. (2.11) holds for this (X, P) and /(f) given 
by Theorem 3.4-. 

Proof. Let (X k , P k ) be the characteristics associated to F k . Their uniform conver- 
gence follows from (2.6-2.7). As the f k are classical solutions to the Vlasov equation, 
they are also mild solutions, i.e. they satisfy 

/ <p(x,p)f k (t,x,p)dxdp = [ ij(X k (t:0,x,p),P k (t- 1 0,x,p)) f (x,p)dxdp (3.18) 
Jr 2 Jm. 2 

for all test functions ip € C C (R 2 ) and all f > 0. 

The sequence ^(X k (t; 0, x,p),P k (t; 0, x,p)) fo(x,p), resp. ip(x,p) f k (t,x,p), con- 
verge uniformly on the support of towards V(X(f; 0, x,p), P(t; 0, x,p)) fo{x,p), 
resp. ip(x,p) f(t,x,p). Hence, the integrals on both sides of (3.18) converge toward 
the two sides of (2.11), and the latter equality is obtained at the limit. □ 

4. Local-in-time existence of characteristic solutions 

In order to have stronger solutions, i.e. such that / satisfies (1.14) in the character- 
istic sense, we need estimates on the Lipschitz constant of the force, which in turn 
amounts to bounding the second derivative of the vector potential. In this Section 
we show that, in both the NR and QR cases, the limiting vector potential given by 
Theorem 3.2 does satisfy such a bound, at least for a short time. 

In the sequel, the operator £ will be that of Subsection 3.1. We no longer consider 
the extended version, so that we can obtain estimates based on characteristics. 
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Lemma 4.1. Let (E, A) e C° (0,T; W^°°{R) x TT 2 ' 00 ^)), and (£',4') = 
£(_E, A). T/ien, ifte second derivative of A' is bounded as: 



d 2 A' 



dx 2 



<C + t 



OA 



dx 



(4.1) 



Proof. We use the formula (A. 4), with / replaced by —nA, and the bounds (2.15) 
and (2.20) for n: 

rt 



d z A' 



dx 2 



(*, x) 



<C + 



<C + t 



OA 



dx 
OA 



\A\ 



dn 



o 
OA 



dx 



\A\ 



dx 
dn 



dx 



(s,x + s — t) ds 
(s,x — s + t) ds 



dx 



(C 1 + C 2 ||^|| t )+t||A|| t C 3 e t ( 1+ l' a ^l' t ). 



□ 



Theorem 4.1. Let (E k , A k ) ke ^ an d C^i-^) ^ e as ^ n Theorem 3.2. There exists 
< T* < T such that, for < t < T* , the sequence (d 2 A k ) k( _ N is uniformly bounded 
in x and t. As a consequence, the couple (E, A) allows to define a characteristic 
solution f to (1.14) on the interval [0,T*). 

Proof. We have: 

d x F k = d x E k + {d x A k f + A k d 2 x A k . 
Thus, the boundedncss results of Theorem 3.2 and (4.1) yield: 

||^F fe+1 || t <a + /3^e*( 1+ II^H0, 

for some positive constants a, (3 depending only on the initial conditions, T and 
(possibly) L. Hence, ||<9 x -Ffe|| t < v k (t), where v k {t) is defined by the recurrence 
formula: 



v (t) = \\d x F a \\ t , v k+1 {t)=a + pte^ 1+Vk ^ := <pt{v k (t)). 



(4.2) 



In Appendix B we show that, for t < T*, the sequence (v k (t)) k <z^ is convergent and 
hence bounded. This gives the uniform boundedness of (d x F k (t)) ke N and, by (4.1), 
that of (d 2 A k (t)) ke fi, on [0, T*) xR. Hence, the latter sequence admits a subsequence 
which converges weakly-* in L°°([0, T*] x R) towards a limit which is necessarily 
equal to d x A. This gives the last part of the conclusion. □ 



5. Global characteristic solutions in the quasi-relativistic case 

In this section, we only consider the QR case. We show that, in this framework, 
characteristic solutions exist for any time. This result rests on a subtler treatment 
of the formula (A. 4). On the other hand, we need to introduce another assumption: 
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Hypothesis 5. The majorising function g of Hypothesis 2 moreover satisfies 

\p\ 2 g(p)dp < oo. 



With this hypothesis, we can derive a better bound on d 2 A. 

Lemma 5.1. Let E E C° (0,T; W 1 '°°(M.)), A e C° (0, T; VF 1,00 (IR)) n 

C 1 (0,T;L°°(R)) ; and (£',A') = £(£,,4). T/ien, the second derivative of A' is 
bounded as: 
8 2 A' 



dx 2 



<Co + dn 3 (*||F|| t ) 





dA 




dA 


{ 


dx 


+ 

S 


dt 



+ ||A|| a (l + ||F|| a ) , (5.1) 



where II3 is a polynomial of third degree, whose coefficients are positive and depend 
only on the function g. 



Proof. We use the formula (A. 4) with / replaced by —nA: 

d 2 A' /"* 
2 ^T^ X ) ^ D (t,x) + j 



dA\ , . / dA\ , , , 

n — — (s,x + s — t) - \n — — (s,x + t- s) ^ ds 



dx 



dx 



dn 



( s > x + s -^ + ^i^J (s,x + t-s)j ds (5.2) 
where D n depends on the initial data only. The first line is bounded by (2.15) as: 



C + (C 1 +C 2 t\\F\\ t ) f 
Jo 



f 


dA 


Jo 


dx 



ds. 



In order to bound the second line in (5.2), we use (1.14) to rewrite 

df d df 

(p-1) g^(s,x + s-t,p) = - — [f(s,x + s-t,p)] + F(s,x + s-t)—(s,x + s-t,p). 

We can integrate this w.r.t. p, since p — 1 never vanishes. Hence, the first part of 
the integral which appears on the second line in (5.2) becomes: 



h 



- f A(s,x + s-t) [ ^— 
Jo JrP-1 



d_ 
"d~s 



[f(s,x + s -t,p)] 



+ F(s,x + s - t) —(s,x + s - t,p) 



dpds ; 



then, performing integration by parts in s and p: 



h = -A{t,x) / — — — dp + A (x-t) / — - — - 
Js. P ~ 1 Jm. P ~ 1 

f* 9 f f(s,x + s- t,p) 

+ / —[A(s,x + s-t)\ / — dpds 

Jo ds J R p-1 



dp 
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+ 



f (AF)(s,x + 8-t) f f(s,x + s-t,p) 
Jo Jr 



dp 



p-l 



When p tends to +oo, there holds: 

1 ~ 2 d r 1 

p — 1 dp 



p-l 



P+^/l+P 2 



^T+p 2 (p - ^T+p 2 ) 



dpds. 



-4p, 



so that the integrations by parts mentioned above are fully justified using Hypothe- 
ses 1, 2, 5 and moreover, 

\h\ <C + Ci||A|| t f (1+ P 2 )f(t,x,p)dp 
Jr 

+ C 2 f {\\d x A\\ s + \\d t A\\ s + \\A\\ s \\F\\ s } f (l+p 2 )f( s ,x + s-t,p)dpd s . 
Jo Jr 

Using Lemma 2.5, we obtain a bound on the integrals in p, so that: 





dA 




dA 


{ 


dx 


+ 

S 


dt 



+ ||A|| fl (l+||F|| 8 ) . 



<Co + CiIl3 (f||F|| t ) 



Of course, L 2 :— J (A^) (s,x + s — t) ds is bounded in the same manner. 
Hence (5.1). □ 

In the case of two solutions corresponding to two pairs of fields (E\,A{) and 
(E 2 , A 2 ), we have similarly: 

Lemma 5.2. Let [E\,A\) and (E 2 ,A 2 ) be as in Lemma 5.1; and [E'^A'A = 
C(E u Ai). We have: 



d 2 A[ 


d 2 A' 2 




8A 1 


dA 2 


+ 


dA 1 


dA 2 


dx 2 


dx 2 




dx 


dx 


t 


dt 


dt 



+ Pi -A 2 \\ t + \\F 1 -F 2 \\ t j (5.3) 
where the constant K is polynomial in (t, \\Ai\\ t , \\d x Ai\\ t , \\dtAi\\ t , \\Fi\\ t ). 
Proof. This time, we have to bound the integral 

as well as a similar one in which the current point is (s,x — s + t), and which will 
be handled in exactly the same manner. Writing: 

djni A x ) _ d(n 2 A 2 ) _ A f _ + ^Hl (A - A) 

dx dx 1 \ dx dx J dx 1 2 



+ ni 



8Ax dA 2 

dx dx 



dA 2 

dx 



(m - n 2 ) 
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we split I into four parts I\, I 2 , A, I4, which are the integrals of the four terms 
above. By (2.15) and (3.10), we obtain: 

|A| < {C + d ||F 2 || t ) t \\d x A 2 - ; 

|A| < (c + d max ||Fj|| t ) \\d x A 2 \\ t t \\F 2 -Fi\\ t . 

Then I 2 is bounded as in the previous Lemma: 

\h\ <Cn 3 (t \\F 2 \\ t ) {H^Ai - dxAzWt + WdtA, - dtAiWt+WA! - A 2 \\ t (1 + ||F 2 || t )} . 

There remains to bound I\. Performing the same computations as in Lemma 5.1, 
we obtain: 



A 



= -A 1 (t,x) [ 
Jr 

■* d 



(h-h)(t,x,p) 



ds 



p-i 

[A-L(s,x + s-t)] 



dp 

(A - f2)(s,X + S-t,p) 
P-I 



dpds 



+ / (A 1 F 1 )(s,x + s-t) f fi(s,x + s-t,p)-^- 
Jo Jr d P 

- [ (A 1 F 2 )(s,x + s-t) [ f 2 (s,x + s-t,p)-^- 
Jo Jr dp 



1 



p-l 
1 



dp As 
dp As. 



We rearrange the last two integrals by writing, as usual, F\ fi — F 2 f 2 = F\ (/i — 
A) + (F 1 -F 2 ) f 2 . Hence: 

\h\<c\\\A x \\ t [ (l+p 2 )(f 1 -f 2 )(t,x,p)dp 
I Jr 

+ I {||a a Ai|| a + 11^11, + ||i4 1 || a ||Fi|| a } / (l+p 2 )(A - f 2 )(s,x + s-t,p)dpi 
Jo Jr 

+ f || A 1 ||a||Fi-F 2 || a / (l+p 2 )f 2 (s,x + s-t,p)ApdsX, 
Jo Jr J 

which is bounded with the help of Lemmas 2.6 and 2.5 as: 

|A| < C |{||0*Ai|| t + \\dtAx\U + \\Ai\U (1 + ||Fi|| t )}Il3 (t max ||^|| t ) J* ||^ - F 2 

+ ||A 1 || t n 3 (t||F 2 || t ) f WF, - F 2 \\ s . 
Jo 



-2\\s 



Putting all these bounds together, we obtain (5.3). 



□ 



Such a bound as (5.3) cannot be derived in the NR case, basically because the 
velocity p — p is not bounded: singularity formation could happen at any speed. 
The proof of the above two Lemmas does not apply, because neither l/(p — 1) nor 
its derivative — l/(p — l) 2 are locally intcgrable near p = 1, and the integrations by 
parts would be unjustified. 
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Nevertheless, in the QR case we immediately deduce the following theorem: 

Theorem 5.1. Let (Ek, Ak)keN an d (E,A) be as in Theorem 3.2. The sequence 
(dlAk) keN converges uniformly in x and t. As a consequence, the triple (f,E,A) 
defines a characteristic solution to (1.1^-1.16) on the interval [0, T}. 

Proof. The convergence of (d 2 Ak) fceN immediately follows from (5.3) and the con- 
vergence results of Theorem 3.2, and the limit is necessarily equal to d 2 A. So, d x Fk 
converges uniformly towards the x-derivative of the force given by (E, A); and we 
get the conclusion. □ 

6. Stability of equilibria 

In this section, we restrict ourselves to classical solutions to the system (1.14-1.19). 
Thus, in the NR case, we take T > small enough to have a classical solution 
on [0,T\. 

Furthermore, we only consider solutions of finite mass, energy and entropy — 
the latter concepts being precised in the following Subsection — in order to rule 
out some unphysical or pathological behaviours. 

6.1. Energies and entropy for the system (1.14 1-19) 

The physical energy of the system naturally splits in two parts, the transversal 
energy made up of the terms containing the vector potential A, and the longitudinal 
energy for the other terms. 

The transversal energy is written WT[f(t), A(t), d t A(t)}, where: 

WT[f,A,A]:= l - J [n[f](x)A(xf + \8 x A(x)\ 2 +A(x) 2 ) dx (6.1) 

and, of course, n[f] denotes the density of /. In the periodic case, the integral is to 
be taken over one space period; in the open-space case, over the whole space. 

Let V = H\R) or #£(K) and H = L 2 {R) or L 2 L (R). The norm denoted || • ||, 
without any subscript, will be that of H; and (• | •) is its scalar product. For any 
function w e C° (0, T; H), we set 

(w) t = sup ||w(s)|| 2 . 

0<s<t 

We have the following result: 

Proposition 6.1. Assume that A G V and A E H. Then A e C°(0,T;V) n 
C 1 (0,T;7i). ,4s a consequence, the transversal energy is finite-valued and differen- 
tiable on [0,T], and there holds: 

±WT[f(t),A(t), 8 t A(t)] = \l d t n(t, x) A(t, x) 2 dx. (6.2) 
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Proof. Consider the approximating sequence (Ak) keN from Theorem 3.2. Clearly, 
A Q G C° (0, T; V) n C 1 (0, T; H). Assuming that A k G C° (0, T; V) n C 1 (0, T; H), we 
see that A^+i can be identified with the variational solution to: 

^(A k+1 (t) \ B) + (d x A k+1 (t) \d x B) = -{n k (t)A k (t) \ B) , VBeV, 
A k+1 (0) = A o e V, d t A k+1 (0) = A o eH, 

where, of course, n k — n[f k ]; recall that, thanks to (2.15), \\n k \\ t is uniformly 
bounded by a constant on the interval [0,T]. Hence A k+1 e C°(0,T;V) n 
C 1 (0,T;7i). By induction, the whole sequence belongs to this space. 
Moreover, we have the classical energy estimate: 



d_ 1 

dt 2 



dA k+1 



dt 



OA 



fe+i 



dx 



n k (t) A k (t) 



dt 



from which we deduce, by Young's inequality: 
1 / dA k+1 \ I dA k+ i 



2 \ dt / t + \ dx 
On the other hand, 



||A fc (t)|| 2 = || J 4 || 2 + 2 



<C + Ci (A k ) t . 



(6.3) 



(6.4) 



o 



A k (t) 



applying once more Young's inequality, we get: 

'dAk 
dt 



{A k {i)) t <C Q + C 1 f 
Jo 



dA k 
dt 



ds. 



ds ; 



(6.5) 



This inequality, together with (6.4), shows that 



+ 



dA k 
dx 



ds, 



t \ — " / t 

which gives the uniform boundedness of (d t Ak) t and (d x A k ) t by Lemma 3.5, and 
that of (A k ) t by (6.5). 

As a consequence, for any t the sequences (A k (t)) keN (d x A k (t)) keN , (d t A k (t)) keN 
admit weakly convergent subsequences in H, hence (A(t) 7 d x A(t) 7 d t A(t)) G H 3 . 
Then, similar computations show that the sequence (^4fc) fegN is indeed strongly 
convergent in C° (0, T; V) n C 1 (0, T; H); and, at the limit (6.3) gives 

d_ 1 

dt 2 

from which we deduce (6.2) by integration by parts thanks to Corollary 2.1. □ 

We now study the longitudinal energy. To this end, we first define the electro- 
static field and potential functionals. In the open-space case, we set 





dA 


2 


dA 




dA\ 


( 


dt 


+ 


dx 


j =-(n(t)A(t) 


~&t) 



E[f] = -dMf] := d x (0 cxt - 0[/]) , 
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where (j) ext is (the opposite of) the external confining potential, satisfying 0" xt = n cxt 
and (j>[f] > is the second primitive of the density n[f] of /, given by the formula: 

mx) = \ I \x-y\n[f](y)dy (6.6) 
z Jr 

due to Lemma 2.8. We remark that is a linear, self-adjoint and positive operator. 
As for the potential 4> ex t, we assume that it satisfies the following hypotheses. 

Hypothesis 6. There exist C > and m e N s.t. 

|0e*t(aO| < C (1 + \x\ m ) and f \p\ m g(p)dp < oo. 

Jr 

(g is the majorising function in Hypothesis 2.) 

Let us remark that, if < n m ; n < n cxt < «max, the above hypotheses are 
satisfied with m = 2, provided Hypothesis 5 holds. 

In the periodic case, we do not use the above decomposition of the potential, 
since the external and the self potential would not be periodic. We set: 

E[f] = -d x *[f\, 'dim = n cx t - n[f], 
the uniqueness of $[/] being ensured by imposing periodicity and 

/ §[f]{x)n elA {x)dx = 0. 
Jo 

This, in turn, ensures the self-adjointness of the operator $. 

As for the kinetic part, we denote by n{p) the primitive of p, viz. 

p 2 , 

k(p) = — (NR case), n(p) = yl + P 2 (QR case). 

We recall that, in the QR case, we always assume Hypothesis 5. In the NR case, the 
following assumption will be needed to ensure the differentiability of the energy. 

Hypothesis 7. The majorising function g satisfies 

/ \p\ 3 g(p)dp < oo. 
Jr 

We are now ready to state the 

Definition 6.1. In the periodic case, the longitudinal energy functional is given 
by the two equivalent formulae: 

WL[f] = £j | K (p) - \ *[/](*)} f(x,p) dpdx, (6.7) 

WL[f]= [ [ K (p)f(x,p)dpdx+l [ |9 x $[/](x)| 2 dx; (6.8) 
Jo Jm. 1 Jo 
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while in the open-space case, one sets: 

WL[f] =11 L(p) +4> ext (x) - l<p[f}(x)\ f(x,p)dpdx. (6.9) 

Let us remark that the definition (6.9) becomes formally equivalent to (6.7) 
or (6.8). However, in the open-space case |i9z0[/]| 2 is clearly not integrable and 
the integration by parts to arrive to a formula like (6.8) are not justified. We will 
comment on this in the next subsection. 

Now, we can analyse the evolution of the longitudinal energy. 

Proposition 6.2. Let (f,E,A) be a classical solution to (1.14-1-19). Assume Hy- 
potheses 1, 2 and 7 if the system is NR, Hypotheses 1, 2 and 5 if it is QR. In 
the open-space setting, assume moreover that Hypotheses 3, 4 an d 6 hold. Then 
1 1 ► WL[f{t)\ is finite-valued and differentiable for all t G [0, T]. Namely: 

±WL[f(t)} = Jj(t,x)^(t,x)dx. (6.10) 

Proof. First, we consider the case of periodic solutions. The fmiteness of the energy 
defined by (6.8) follows from Lemma 2.5 and the various boundedness results of 
Theorem 3.2. By the way, the boundedness of E yields that of $; and integrating 
by parts gives the equivalent expression (6.7). Then, as in Corollary 2.1, we bound 

d df df 

— [ K (p)f(t,x,p)] = -p K (p)^L + K (p)F(t,x)^-, 



d 

— [K(p)f(t,X,p)} 



<k{ p ) (F, + \p\)gti\ F \u(p)^ {1+l ^ Fh) - 



The highest power of p is p 3 (NR) or p (QR): the above function is integrable 
in p, and in x given the finite length of (0, L). Finally, Eq. (1.15) gives d t {\ E 2 ) = 
Ej: these functions are continuous on (0,L), so the second term in (6.8) is also 
differentiable. 

Let us now consider the open-space setting. The fmiteness and differentiablity 
of the integral JJ {n{p) + (f> cx t(x)} f(t,x,p)dxdp follows from Lemma 2.9. Then, 
the self potential defined by (6.6) satisfies: 



0<4f(i)](z) < \ [\x\ j n[f](y)dy + j \y\n[f](y)Ay 



Z being defined in Lemma 2.8. By the same token, 

0< JJ <j>[f(t)}(x)f(t,x,p)dxdp< \ J(M\x\ + Z)f(t,x,p)dxdp<MZ. (6.12) 

There remains to check the differentiability of the self potential energy Using the 
Vlasov equation (1.14) and Corollary 2.1, we find: 

\d t f(t, x,p)\< (F* + \p\)g t || F |, (p) g t2 (x - tp) e 4 (i+II^^IU). (6 .1 3) 
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But the linearity and sell-adjointness of <j> imply: d t ((j>[f] f) = 2d t f<fi[f], which 
allows to bound 

j^[/(*)]/(*)](s,P) < (Co + CilxQ (l + F* + \p\)g tllFh (p)g t2]lFh (x-t P ) 



by (6.11) and (6.13). This proves the integrability of this function, and the differ- 
entiability of the self potential energy. 

Finally, Eq. (6.10) is obtained through tedious but straightforward computations 
(cf. Propositions 1.5 and 1.6 in 2 ), all the integrations by parts being justified by 
the above arguments. □ 

From (6.2), (6.10) and the continuity equation (2.19), we immediately deduce: 

Theorem 6.1. Under the hypotheses of Propositions 6.1 and 6.2, the energy 
W(t) := WT[f(t),A(t),d t A(t)]+WL[f(t)] is constant. 

For reference, we give the formulae for the energy in the FR case. There is no 
splitting in transversal and longitudinal parts. The total energy is given by: 

<- L r r i i 

f(t,x,p)dpdx 



/o 

+ 



w{t) -=U 

2 Jo 



y/l+p* + A(t,x)*--$[f(mx) 



dA 
~dt 



(*, x) 



dA 

dx 



(*, x) 



dx 



in the periodic case, and 

W(t):= J^j y/1 +p 2 + A(t,x)i + Mx) - \ 



+ 



z Jr 



dA 
~dt 



(t,x) 



+ 



dA 

dx 



{t,x) 



dx 



(6.14) 

f(t,x,p)dpdx 

(6.15) 



in the open-space case. The reader may check that these expressions are formally 
constant, even though the mere existence of classical solutions is an open problem, 
and subtler arguments are probably necessary to justify the calculations. 

Definition 6.2. Let a £ C 2 ((0,+oo)) n C°([0, +oo)) be a strictly convex and 
bounded-from-below function, which satisfies 

a(s) 



lim 



+00. 



s^+oc s 

Let 7 denote the generalised inverse of — a' (extended by if necessary): it is a 
decreasing function in its support. 

The entropy associated to a distribution function / is defined as: 



8,1/]-//. 



(f(x, p)) dp dx. 



The most classical example is a(s) = s Ins — s, i.e. j(s) — e s , associated to 
Maxwellian distribution functions. 
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Clearly, S a is a convex, bounded-from-below, weakly lower semicontinuous func- 
tional on its domain of definition. From 5 , we know that for any /* e _L 1 (R d ), there 
exists a function a as above s.t. is finite. Thus, we shall choose a according 

to the following 

Hypothesis 8. In the periodic case, J a(g(p))dp < oo. In the open-space case, 
J o-(g(0) g(p)) dp < oo and JJ a(g(x) g(p)) dx dp < oo. 

In this case, the following identities hold: 

J a {g r {p)) dp = 2ra(g(0)) + J a(g{p))dp; 

J J a(g p (x)g r (p)) dxdp = 4prg(0) + 2(p + r) J a{g(0) g{s)) ds 

o(g{x)g{p))dxdp. 



Then, the finiteness and differentiability of S a [f(t)] stem from arguments very 
similar to the proof of Proposition 6.2, and one can easily deduce. 

Theorem 6.2. Under Hypothesis 8, the function S a [f{t)\ is constant on [0,T]. 
6.2. Equilibria of (1.14-1.19) 

We are looking for solutions to (1.14-1.19) which do not depend on time, i.e., 
solutions to the coupled problem in (f(x,p),A(x)): 



p^ + — 

dx dx 



A 2 ' 



-^ + n[f]A = (6 ' 17) 

with the potential satisfying the Poisson equation — -j^z<&[f] = n cxt — n[f\. 
We have the following result 

Lemma 6.1. Any classical solution to (6.16-6.17) s.t. M > and A £ V is a 

Vlasov-Poisson equilibrium, i.e. it has the form (/oo,0), where f^ solves 

^ + ^1^=0, (6.18) 
ox dx up 

with — -j^z<&[f oo\ — n cv ± — n [foo\ and consequently satisfies: 

/oo(*,p) = -F(«(p) - *[/«,](*)) , (6.19) 

for some function T . 

Proof. If A G V, then for any B eV, 

a f (A, B) := j (A'(x) B'{x) + n[f]{x) A{x) B{x)) dx = 0. (6.20) 
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Now, take B = A to deduce that a,f(A, A) = 0. Since both terms in a,f(A, A) are 
non-negative, we deduce that 

J \A'{x)\ 2 dx = j n[f](x) A{x) 2 dx = 0. 

From the first identity we conclude A is constant, while from the second we conclude 
A = since n[f] > on some non- negligible subset of R or (0, L). We are left 
with (6.18), whose solution is well-known to be of the form (6.19). □ 

The function T can be precised by demanding that the solution should minimise 
a "free energy" functional. In other words, the choice of the entropy function a is 
determined by the particular equilibrium one is interested in. 

Definition 6.3. Let / e ^(R 2 ), rcsp. / G L^(M 2 ), the free energy of / is: 

K„[f\ :=WL[f] + S a [f]. 

Let then A e V and A G H; we set 

KT a [f,A,A] :=K a [f]+WT[f,A,A}. 

We consider the set of suitable distribution functions with fixed mass, i.e. 

K{L,M) := {/ £ L\(R 2 ) : / > a.e. and ||/|| L i (R 2) = m} . 

Lemma 6.2. In the periodic case, KT a is a strictly convex and bounded- from-below 
functional on IC(L, M) x V x H. It has a unique global minimum which takes the 
form (/oo,cr,0, 0), where 

fco,a = 7(«(P) - $[/«,,„] (*) - ") (6.21) 

and is therefore a stationary solution of the Vlasov-Poisson system. The constant 
a is uniquely determined by M and L. 

Proof. Eq. (6.7) shows that WL is convex in this case; it is clearly lower semicon- 
tinuous and on /C(L, M) it is bounded from below. Since S a enjoys the same prop- 
erty and is strictly convex, we deduce that K a has a unique global minimum ./over- 
writing the Lagrange equation expressing the minimisation under the constraint 
J f = M (cf. 3 ' 5 ), yields the formula (6.21), where a is the Lagrange multiplier. 
Then, it is clear that for any (/, A, A) ^ (foo,a, 0, 0), 

KT„[f, A, A] > KT a [/, 0, 0] > KT a [f^, 0, 0] , 

where the last inequality is a consequence of the results in 3 ' 5 . □ 

A similar result in the open-space case does not hold in our case, in contrast 
to the situation studied for nonlinear stability of the Vlasov-Poisson system in 
higher dimensions. 17 ' 4 ' 3 < 5 The main difference in ID being that d x (f>[f] £ L 2 (R) as 
pointed out before. Sobolev embeddings in d > 2 allow, under confining conditions 
on the external potential e xt 7 to deduce a result similar to previous Lemma. 10 ' 3 ' 5 
Moreover, the functional KT a ceases to be convex in our ID open-space case. 
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6.3. L p -nonlinear stability of equilibria in the periodic case 

In this subsection, we just collect the known results in several references and applied 
to the particular case we deal with. Like in 17 - 4 ' 5 , one can rewrite: 



KT a [f, A, A] - KT a 0, 0] = S CT [flf^] + WT[f, A, A], (6.22) 



where the relative entropy of the distribution function / w.r.t. g is defined as: 



Yet, as a consequence of Theorems 6.1 and 6.2, KT a [f(t),A(t),dtA(t)] is con- 
stant for any classical solution to (1.14-1.19). This implies (see 17 ' 4 ' 5 ) that: 

(1) The LP norm of / - is bounded for 1 < p < 2, if infer" {s)/ S p- 2 > 0. 

(2) The L 2 norm of /— /oo,o- is bounded, if ./oo,o- is a Maxwellian, i.e. a(s) = s In s—s. 

(3) The H 1 norm of $[/] - is bounded. 

(4) The transversal energy is bounded — indeed, we already knew this, and even 
a little bit more, from Proposition 6.1. 

Here, all the norms are taken for x £ [0, L], p £ M. 

The first three points follow from §§3 and 4 of 5 and references therein. Indeed, 
the arguments in these passages are independent of the dimension. 

In other words, the one-dimensional periodic Vlasov-Poisson equilibria are L p - 
nonlinearly stable under one-dimensional Vlasov-Maxwell perturbations. 

Let us finally mention that Landau damping was proved in 6 in the case of the 
Vlasov-Poisson system in the periodic case. As a consequence, it was proved in 6 
that some Vlasov-Poisson equilibria which are L 1 -nonlinearly stable, are unstable 
for a weak topology. This is not known to happen in our ID Vlasov-Maxwell system, 
although numerical computations seem to indicate that nonlinear Landau damping 
should occur in this model. Let us point out that there is no contradiction between 
these two stability assertions, since weak topology neighbourhoods of the equilibria 
are much larger than L 1 neighbourhoods. 

Appendix A. The Duhamel formulae 

The following representation formulae are capital in the various estimations of Sec- 
tions 3-5. 

The unique temperate solution u to the wave equation 




I u(0,a;) = u (x) £ X, d t u(0,x) = v (x) £ X, 
(K) n <S'(R), is explicitly given by the formula: 




where X = L 



loc 
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+ / / f( Sl y)dyds\. 

JO Jx+s-t ) 



(A.l) 



Hence, after some computations, the formulae for the various derivatives of u: 



du . 



dt 



du 



1 



-^r{^ x ) = o { u o( x + *) - u 'o( x ~ *) + v o( x +t) + v (x - i) 



+ J [f(s,x + s-t) + f{s,x + t-s)]d s y (A.2) 



— (t, x) = - yi' {x + t) + u' Q (x - t) + v (x + 1) - v (x - t) 
+ J [-f{s,x + s-t) + f{s,x + t-s)}ds^, 

d 2 u 1 f 

g^(t,x) = -iv%(x + t)+v£(x-t)+r/ (x + t)-v' (x-t) 



+ 



df , ^ df , 

~z-(s,x + s- t) + —{s,x + t- s) 
ox ox 



ds 



which are valid e.g. if u' , u' \ v' G X and / e C 1 (0, T; X). 



(A.3) 



(A.4) 



Appendix B. Convergence of the sequence defined by (4.2) 

Here is the technical lemma announced in the proof of Theorem 4.1. 
Lemma Appendix B.l. Let {vk{t))km an d ft be defined by (4-2). 

(1) There exists T\ < +oo such that, for < t < T\, ip t admits two fixed points 
v t <v*. Ifvo(t) < v l , the unstable fixed point, then the sequence Vk(t) converges 
toward the stable point Vt ■ 

(2) For t > Ti, there is no fixed point, and (vk(t))kew diverges to +oo. 

(3) When t tends to zero, v l goes to infinity, while v t remains bounded. 

Consequently, there exists < T* < T\ s.t. (vk(t))keN is convergent for < t < T* . 

Proof. As (fit is convex and increasing, there are only three possibilities: 

(1) The equation <ft{x) = x admits two solutions v t < v*, which are respectively 
stable and unstable fixed points: < (p' t (v t ) < 1 and <^(v*) > 1. 

(2) The two solutions merge in a unique fixed point v, which satifies both (ft(v) = v 
and (fi' t (v) = 1; it is stable on the left side, unstable on the right side. 

(3) There is no fixed point at all, and ifit(x) > x for all ieR. 

The case 1 is achieved for t small enough. Indeed: 

lim/3ie t(1+a:) = 0, VxeM+. 
t->o 

Hence, for any fj, > 1, there exists > s.t. ip t ((j,a) < /ia for t < t^. On the other 
hand, <Pt{o() > a. Hence, ip t has (at least) one fixed point in the interval (a,/ia) 
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when t < tp. Given that <fit{%) 3> x when x — > +oo, there is another fixed point in 
(fia, +oo). 

Then, we notice that, when a; is fixed, <£>*(x) is a decreasing function of t. This 
has two consequences. Firstly, if s < t and ipt has fixed points, then ip s also has 
fixed points. This proves the existence of Ti, which may be finite or not. Clearly, if 
Ti is finite, it achieves the case 2. Eliminating v between (ft{v) — v and <p' t (v) = 1 
gives the following equation for t: 

(3t 2 e at2+2t = 1. 

As the l.h.s. is zero for t = 0, infinite for i = +oo and strictly increasing in t, the 
equation admits a unique solution T\ G (0, +oo). 

The behaviour of Vk+i(t) = <fit(vk(t)) then follows from the elementary theory 
of sequences (Figure 1), and the claims 1 and 2 are obtained. 




Fig. 1. Dynamic of the sequence Vk+i(t) = ipt (ffc(*))- 



The second consequence is that v t and v* are resp. decreasing and increasing 
functions of t. Hence, Vt is bounded when t — > 0. On the other hand, i>* > x*, where 
x t is defined by f' t { xt ) — 1. The latter equation gives 

/3t 2 e t(1+a;t) = 1 x* = - t In ^-^^ - 1 -> +oo when t 0. 

This proves the third claim. Finally, let us remark that vq {t) is indeed independent 
of t: E and A , hence F are constant in time. Consequently, for t small enough, 
v f > vo{t). This gives the last part of the conclusion. □ 
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